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ABSTRACT 

In the present work we propose an innovative estimation method for the minimum 
Doppler factor and energy content of the 7-ray emitting region of quasar 3C 279, using 
a standard proton synchrotron blazar model and the principles of automatic photon 
quenching. The latter becomes relevant for high enough magnetic fields and results 
in spontaneous annihilation of 7-rays. The absorbed energy is then redistributed into 
electron-positron pairs and soft radiation. We show that as quenching sets an upper 
value for the source rest-frame 7-ray luminosity, one has, by neccessity, to resort to 
Doppler factors that lie above a certain value in order to explain the TeV observations. 
The existence of this lower limit for the Doppler factor has also implications on the 
energetics of the emitting region. In this aspect, the proposed method can be regarded 
as an extension of the widely used one for estimating the equipartition magnetic held 
using radio observations. In our case, the leptonic synchrotron component is replaced 
by the proton synchrotron emission and the radio by the VHE 7-ray observations. We 
show specifically that one can model the TeV observations by using parameter values 
that minimize both the energy density and the jet power at the cost of high-values of 
the Doppler factor. On the other hand, the modelling can also be done by using the 
minimum possible Doppler factor; this, however, leads to a particle dominated region 
and high jet power for a wide range of magnetic field values. Despite the fact that we 
have focused on the case of 3C 279, our analysis can be of relevance to all TeV blazars 
favoring hadronic modelling that have, moreover, simultaneous X-ray observations. 
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1 INTRODUCTION 

Blazars, a subclass of Active Galactic Nuclei, emit non- 
thermal, highly variable radiation across the whole electro- 
magnetic spectrum. According to the standard scenario, par- 
ticles accelerate to relativistic energies in the jets of these 
objects which point, within a small angle, towards our di- 
rection and the resulting photon emission is boosted due to 
relativistic beaming. 

Detailed modelling of observations, especially in the 7- 
and X-ray regimes, makes possible the estimation of the 
physical parameters of the emitting region. Thus quantities 
like the source size, the magnetic field strength, the bulk 
Lorentz factor and the particle energy density can nowa- 
days be routinely calculated. Furthermore, from the values 
of these quantities one could obtain meaningful estimates 
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for the particle and Poynting fluxes and use them to make 
comparisons, for example, to the Eddington luminosity of 
the source, connecting thus the black hole energetics with 
the jet power. 

One major uncertainty of the modelling is the na- 
ture of the radiating particles. While there seems to 
be a consensus that the emission from radio to X-rays 
comes from the synchrotron radiation of a population 
of relativistic electrons, there are still open questions 
regarding the 7-ray emission of these objects. Broadly 
speaking, the models fall into two categories, the lep - 
tonic ones (e.g. iDermer. Schlickeiser. fc Mastichiadi s! Il99 
Maraschi. Ghisellini. fc Celottilll992l ; lDermer fc Schlickeise: 
1993ft which assume that the same electrons which radi- 
ate at lower frequencies via synchrotron produce the 7- 
rays by inverse Compton scatt e ring and the hadronic ones 
llMannheim fc Biermam] Il992l ; IMiicke fc Protheroel 12001 : 
iBottcher. Reimer. fc Marschej 120091 : IMiicke et all l2003h 
which postulate that an extra population of relativistic pro- 
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tons produce the high energy radiation as a result of hadron- 
ically induced and electromagnetic processes. 

Due to the very different radiating mechanisms in- 
volved, the two classes of models can result in very dif- 
ferent parameters for the source. Thus, for example, typ- 
ical leptonic synchrotron self-compton models require low 
magnetic field strengths, ranging fr om B ~ 0.01 — 1 G 
for high synchrotro n-peaked BL Lacs (jTavecchio et al.ll201ll : 
iMurase et alll2012( ) up to B ~ 1 - 10 G for Flat Spectrum 
Radio Quasars, and low jet power (~ 10 47 erg/sec) (e.g. 
ICelotti fc GhisellinilEioih . On the other hand, the corre- 
sponding values for the hadroni c models are higher by at 
least one order of ma gnitude (e.g. lProtheroe & Miickc 2001; 
iBottcher et al.ll2009T ). While both can fit reasonably well the 
observations, they both face some problems. For instance, it 
has been argued that one problem with leptonic models is 
the high ratio between the required relativistic electron en- 
ergy density to the magnetic one, implying large departures 
from equipartition. On the other hand, the hadronic models 
imply jet powers which can be uncomfortably high, espe- 
cially when compared to the accretion luminosity. 

In previous work (|Petropoulou fc M astichiadis 20 1^) we 
have argued that automatic 7-ray quenching, i.e. a loop of 
processes which result in spontaneous 7— ray absorption ac- 
companied by production of electron-positron pairs and soft 
radiation, can become instrumental in the modelling of high 
energy sources. Its application to any 7-ray emitting region, 
makes it relevant to both leptonic and hadronic models; the 
fact, however, that quenching requires rather high magnetic 
fields makes it more relevant for the latter. In the present 
paper we revisit the hadronic model taking into account 
the effects of non-linear photon quenching. This, as we shall 
show, has as a result to exclude many sets of parameters, 
which otherwise would give good fits to observations, as the 
non-linear cascade growth modifies drastically the produced 
multiwavelength spectrum. 

As a typical example we focus on the quasar 3C 279. 
This has been detected for the first time in very high- 
energy (VFffi]_j^ys_at > 100 GeV by the MAGIC tele- 
scope l)Albert et al.ll2008l ) and ever since it remains the most 
distant VHE 7-ray source with a well measured redshift. 
Both leptonic and (lep to) hadronic models have been applied 
l|Bottcher et al.ll2009T ) and it has been pointed out that the 
former require the system to be well out of equipartition. 

In the present work we fit only the VHE part of the spec- 
trum using a proton-synchrotron blazar model and use th e 
contemporary X-ray data (see e.g. IChatteriee et all (|2008)) 
as an upper limit; the modelling of the complete multiwave- 
length spectrum requires a primary leptonic population radi- 
ating in the IR/X-ray energy bands which, we assume, that 
can always be determined. We show that one can derive a 
lower limit for the Doppler factor of the radiating blob by 
just combining (i) the values of the proton luminosity and 
the Doppler factor for a given magnetic field that provide a 
good fit to the MAGIC data and (ii) the fact that for high 
values of the proton injection compactness the absorption 
of 7-rays becomes non-linear leading to an overproduction 
of softer photons and consequently violating the X-ray ob- 
servations even in the absence of a leptonic component. We 
also show that the minimum possible value of the Doppler 
factor is largely independent of the magnetic field and this 



has some interesting implications for the energetics of 3C 
279. 

The paper is structured as follows. In §2 we give some 
simple, qualitative estimates on the effects of quenching on 
fiting the multiwavelength data of blazars, in §3 we apply 
these, using a numerical code, to the case of 3C 279 and 
we conclude in §4 with a discussion of our results. For nu- 
merical results we adopt working in the ACDM cosmology 
with Ho = 70 km s" 1 Mpc" 1 , Q, m = 0.3 and Ov = 0.7. The 
redshift of 3C 279 z = 0.536 corresponds to a luminosity 
distance Dl = 3.08 Gpc. 



2 ANALYTICAL ESIMATES BASED ON 
AUTOMATIC QUENCHING 

In this section we will justify the existence of a minimum 
Doppler factor using analytical expressions in the simplest 
possible framework. The functional dependence of various 
physical quantities on the magnetic field strength, such as 
the particle density and the observed jet power, will also be 
derived. 

We consider a spherical blob of radius R moving with 
a Doppler factor S with respect to us and containing a mag- 
netic field of strength B. We further assume that ultra- 
relativistic protons with a power law distribution of index 
s are being constantly injected into the source with a rate 
given by 



?p = Qpo7p ^(yp - 7min)#(7n 



7pJ 



(X) 



where 7 m i n and 7 max are the lower and upper limits of the 
injected distribution respectively and H(x) is the Heaviside 
function. Q po is the normalization constant and is also di- 
rectly related to the proton injection compactness as: 
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where t CI = R/c. This can be further related to the proton 
injected luminosity by the relation 



(4) 



where err is the Thomson cross section. 

Protons will lose energy by synchrotron radiation, pho- 
topair and photopion proce s ses a s has been described in 
iPetropoulou fc Mastichiadis (2012|) and therefore the pro- 
ton distribution function is given by the solution of a time- 
dependent kinetic equation that includes particle injection 
in the form of equation {TJ in addition to particle losses 
and escape. Furthermore, since the loss processes will cre- 
ate photons and electrons, one has to follow the evolution of 
these two species by writing two additional kinetic equations 
for them. The solution of the system of the resulting three 
coupled partial integrodifferential equations gives the corre- 
sponding particle distribution functions and the multiwave- 
length photon spectra can be calculated in a straightforward 
manner. The picture above is rather complicated and it can 
be treated only numerically (see section 3). 

In the case, however, where the cooling of protons does 
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not play a significant role in comparison to particle escape, 
the steady-state proton number density is then simply given 
by 



n p(7) = *crQpo7 S >f° r 7min < 7 < 7n 



(5) 



where we have assumed that the proton escape timescale 
equals the crossing timescale. Moreover, since we fit only 
the VHE 7-rays and not the whole multiwavelength spec- 
trum, there is no need in calculating the evolution of the sec- 
ondary leptons that will radiate at lower energies. MAGIC 
data lie above 10 25 Hz. The 7-ray spectrum in vF v peaks at 



100 GeV or equivalently at frequency v° 



2 x 10 2 



which we will try to fit with the proton synchrotron 
spectrum emitted by the particle distribution of equation 
(JSJ). We will focus on the case where the distribution of 
protons is not cooled and the synchrotron spectrum peaks 
at the maximun synchrotron energy, which implies that 
i C ooi(7max) < t C T- Thus, our analysis is valid under the pre- 
requisite 



-B 2 7 max < 10 3 R in cgs. 

In this framework we obtain our first relation 



6° bs = dSB-fi 
where e° b 



(6) 



(7) 



synchotron power per frequency in the blob frame for the 
proton distribution of equation ([5]) is given by 



Pu = C 2 N P B — 
where 

2e 7 / 2 



Ch 



67rm^ 2 c 5 / 2 



(8) 



(9) 



and N p is the total number of radiating protons in the 
blob,_which is given by N p = icrV<2po(l/7min - l/7max) ~ 
tcrVQ po , where the approximation 7 m ! ix << 1 has been 
made; here we have used also 7mm = 1 and s = 2. This 
can also be expressed in terms of proton energy density u p 
as 



N = =PI . 

P ln(7max)m p c 2 



(10) 



In the above expressions V w irR is the volume of the 
spherical blob. 



2.1 Derivation of the minimum Doppler factor 

By combining equation (J7J) with the integrated power over 
all synchrotron frequencies (equation[8]) one obtains the total 
comoving synchrotron power 



L = 2C 2 N p (v° hs ) 1/2 B 3/2 5- 1/2 



(11) 



that is related to the total observed 7-ray luminosity _L° bs « 
10 48 ers/s by the usual relation L° bs = S 4 L. This leads to 
our second relation 



L; bs = C' 2 u p B 3/2 S 7/2 , 



(12) 



1 Throughout the present work quantities with the index 'obs' 
will refer to the observer's frame, whereas all other quantities 
refer to the blob frame. 



where C 2 = 2C 2 V(v° bs ) 1/2 / ln(7 ma x)m p c 2 . For a specific 
magnetic field we find a relation of proportionality 

2/7 , (13) 



jc -2/7 

OLUt, oc 



which is also verified by the detailed numerical treatment 
(see section 3.2). We note that the proportionality at the 
right hand side of relation (|13[) holds only if the magnetic 
field is not strong enough to cause significant proton cooling 
due to synchrotron radiation. 

The maximum proton Lorentz factor has remained so 
far undetermined. However, the fact that the proton gyro- 
radius r E sho uld be less or equal than the size of the blob R 
(Hillas 1984) sets a strong upper limit 



eBR 



(14) 



where k is a scaling factor that takes values less or equal to 
unity. It has been introduced in order to take into account 
the fact, that good fits can be obtained for certain param- 
eter sets with 7max << eBR/m p c 2 ; as an example see the 
parameters used for figure [4] By combining the above rela- 
tion with equation (J7J we find a lower limit for the Doppler 
factor 



hv° hs and Ci = he/m p c(l + z). The total 8 > 8mm = C 3 B~ 



where 



C* 3 = 



Ci 



2 obs 



R 2 



(15) 



(16) 



Automatic quenching of 7-rays was irrelevant to the 
derivation of the relations presented up to this point. Actu- 
ally, it will not play any role in the evolution of the system if 
the magnetic field is weak enough. There is, in other words, 
a necessary but not sufficient condition for the operation of 
automatic quenching, the so-called feedback criterion. This 
can be derived from the requirement that the magnetic field 
is strong enough so that the synchotron photons of the auto- 
matically produced pairs lie above the thre shold for further 
photon-photon absorption on the 7-rays ([Stawarz fc Kirkl 
120071 ; iPetropoulou fc Mastichiadisll201ll . [2012h . This require- 
ment sets a lower limit for the magnetic field 



B > 8B C 



(17) 



where B cr = 4.4 x 10 13 G. If we set e 7 



\z + 1)6- 



and use equations Q and (|14[1 we find an expression for 
the magnetic field B q , below which the feedback criterion 
is not satisfied. This involves only physical constants except 
for the size of the blob R and the scaling factor k: 



1/10 / m p m c c 



lie 



3/10 



enR 



3/5 



(18) 



In order to have an estimate, the above expression gives 
B q w 3.5 G for R = 3 x 10 16 cm and k = 1. Note that for 
a given R this is also the minimum value of B q . One could 
also express B q in terms of 7max as 



B a 



a/4 



3 \ 3/ 4 

m p m c c J \ _ 3/2 
he ) 7max ' 



(19) 



as long as 7 max does not exceed the value given by equation 
()14|) for k — 1. From equation (|12[) becomes evident that if 
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the Doppler factor decreases then the proton energy density 
or the proton compactness equivalently should increase in 
order to obtain the same observed luminosity. If, however, 
the feedback criterion is satisfied, then the proton compact- 
ness is bounded from above. This immediately sets a lower 
limit for the Doppler factor. We proceed next to obtain an 
expression of <5 m i n . In previous work we have obtained an 
analytical expre ssion for the critical 7-ray compact ness (see 
equation (34) of lPetropoulou fc Ma stichiadis] |201ll )l. which 
will prove very useful for our analytical calculations 



1/2.-1/2 



(2e° hs /a 2 B c 



lV2 



(20) 

and o"o = 4/3 is a nor- 



where C4 

malization constant. The 7-ray compactness is defined as 
l 1 = Lctt / iir Rm e c 3 . Using equation expressed in terms 
of proton compactness £ p instead of total proton number jV p 
we can rewrite the above definition as 



c 5 e p B 3/2 5- 



1/2 



where 



C 5 = 



&C2V{v° hs ) 1/2 
4:TTR 2 m c c 3 ln(7 max ) ' 



(21) 



(22) 



In order to avoid an overproduction of soft photons that 
would violate the upper limit set by the X-ray observations 
we require that £ 7 ^ A^c" J > where A is a numerical factor 
between 1 and 10; the exact value can be estimated only 
numerically. This requirement defines a maximum proton 
compactness 

c t 



x c; B - 



(23) 



which, if inserted in equation (|12p . gives the minimum 
Doppler factor 



(7 T 



AirR\rn c c 3 C4 



■2/T 



B 



-1/7 



(24) 



We note that <5 m i n given by the equation above is a rather 
robust limit, since it depends weakly on the physical quanti- 
ties. We note also that the numerical factor A appears in the 
above expression raised to the power —2/7 and therefore it 
does not affect the value of 5 m i n severely. 

Summarizing, we have found that for B < B q , the min- 
imum value of the Doppler factor is defined by the require- 
ment that the gyroradius is less than the blob size and it 
has a very strong dependence on the magnetic field, i.e. 
<5 m in oc B~ 3 . On the other hand, for B ~£ B q , it is defined by 
the requirement that the proton compactness cannot exceed 
a critical value. The dependence on the magnetic field is in 
this case weak, i.e. S mln oc B~ x ^' ' . 



2.2 Conditions for equipartition 

Estimations of the equipartition magnetic field of a rela- 
tivistically moving blob are common in the literature and 
are based on fitting the low-frequency synchrotron spectrum 
of a power-law dis tribution electrons to radio observations 
l|Pacholczvklll970l ). Here we derive an expression for the 
equipartition magnetic field in a proton-synchrotron blazar 
model, replacing the power-law distribution of electrons with 
protons and the radio with VHE 7-ray observations. From 




1.2 1.4 1.6 



Figure 1. Minimum Doppler factor as a function of B (dashed 
lines) and locus of points satisfying the equipartition condition 
(solid line) for L° bs = 10 48 erg/s, u° hs = 2x 10 25 Hz, R = 3x 10 16 
cm and A = 5. 



equation (|12[) we find that the comoving proton energy den- 
sity is given by 

-3/2^-7/2 



AB~ 



(25) 



where A — L° B C2 ■ Assuming that the energy density 
of secondary leptons is negligible with respect to that of 
protons, we can compute the equipartition field _B oq 

B cq = (STvAf^S- 1 



(26) 



718 L%4&R 



i 6 6/7 K b 2 S 5)" 1/7 ln(7 m -,io)5- 1 inG. (27) 



We note that here and in what follows the convention 
Qx = Q/W x in cgs units was adopted unless stated other- 
wise. The above expression apart from the numerical con- 
stant, is identical to that found using equipartition argu- 
ments i n the context of a lept o nic m odel (see e.g. equation 
(A7) of H arris fc Krawczvnskl <|2002l) ). In the previous sub- 
section it was found that the Doppler factor of the blob has a 
minimum value <5 m i n and its functional dependence on B was 
also derived. One can therefore investigate whether or not 
the system can achieve equipartition under the requirement 
that S = <5 m i n and in this case estimate B eq using equation 

urn 

Figure [1] shows <5 m in(-B) given by equations (|15|l and 



(f24)l and the line 5 oc B cq for R 



3 x 10 ib cm and A = 5. 



The discontinuity occurs at B q , where the functional depen- 
dence of (5 m i n on B changes. The expressions derived should 
not be trusted in the neighbourhood of B q . The points of 
intersection correspond to the equipartition magnetic field 
values. The physical system can be found in equipartition 
either moving with a high Doppler factor and being weakly 
magnetized or moving with a modest Doppler factor and 
containing a stronger magnetic field. 

A robust result of our treatment, which also differenti- 
ates it from other works, is that for a given magnetic field 
there is a minimum value for the Doppler factor of the flow 
that is set either by the gyroradius or by automatic quench- 
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ing arguments. Whenever the analysis of a physical problem 
leads in the derivation of an extremum for some parameter, 
it is interesting to study the properties of the physical sys- 
tem in this case. For this reason, in the following we present 
how the particle energy density and the power of the jet de- 
pend on the magnetic field in the particular case of 8 — <5 m m- 
Using equations (|15|l , (|24|l and (125[) we find that 

u p ((5 min (S)) = w* oc B 9 for B < B q (28) 

and 

Mp(5min(B)) = Wp oc B" 1 for B > B q . (29) 

Note the large change in the dependence of the proton en- 
ergy density on B. The last proportionality is not valid for 
very high values of the magnetic field, since in this case the 
synchrotron cooling becomes important and is against our 
initial assumptions. 

2.3 Jet power estimates 

In an one-zone homogeneous model the emitting plasma is 
confined in a blob with radius R moving with a velocity 
j3c and bulk Lorentz factor F at an angle 8 with respect 
to the line of sight. The radiation is further assumed to 
originate from a region in the jet of volume V ftt ttR 3 . The 
energy densities of particles, magnetic field and radiation 
are contributing to the jet power. Using the approximations 
8 « 1/r, 8 w r and /3 « 1 the observed jet power is given 
by 

-Pjot S = ^R 2 8 2 /3c(u B + Up ar t + Urad), (30) 

where u palt ~ u p , f3 — \/l — 1/r 2 and ii ra d is the radia- 
tion energy density, which for simplifying reasons will not 
be included in the following analytical calculations. It will 
be however taken into account in our numerical treatment 
of section 3. Note also that in our analytical formulation 
we have not included any population of 'cold' protons and 
electrons, which are physically implied by charge conserva- 
tion arguments. However, their c ontributions to the total jet 
power are minor (see for example lProtheroe fe Muckell200lh 
and would not affect our results significantly. 

Inserting the expression for u p given by equation (|25|) 
in the above relation we find that the jet power is a function 
only of the product 5 ■ B 

P£ s = kR 2 c (&2f- + A{6B)~*A . (31) 

The jet power is minimized when the following condition 
holds 

B = (6tvA) 2/7 T 1 . (32) 

Comparison of equations ()26[) and (|32[) shows that whenever 
the system achieves equipartition the jet power also mini- 
mizes at the value 

<U, = 5.8x 10 47 « 4 S 8 ln( 7max , 10 ))? f ] erg/ 8 - ( 33 ) 

Thus, even the minimum jet power obtained in the frame- 
work of the proton-synchrotron blazar model is high com- 
pared to values inferred from pure leptonic modelling of 



blazar emission; see e.g. ICelotti fc Ghisellinil (|2008l ) . Since 
for a large range of B values the system lies far from equipar- 
tition (see figure [l|, equations f|2T|) and (|32[) imply that the 
calculated jet power differs significantly from its minimum 
value. 

The fact that the Doppler factor has a lower limit im- 
plies that for magnetic fields higher than a certain value B^ , 
equation (|32|l cannot be further satisfied and there is no set 
of parameters that lead to minimization of the jet power. 
The characteristic value B^ can be derived by combining 
equations (J24J and (|32|l . It is given by: 

B b = 11.6 7?r 6 2/3 A 1/3 ln( 7max , 10 ) 1/3 G. (34) 

The same holds also for weak enough magnetic fields. Thus, 
by combination of equations (|15fl and (|32[) one derives a 
second characteristic value B\ for the magnetic field, below 
which there is no set of parameters that could minimize the 
power of the jet 

Bi = 0.96^6 4/7 (^2 S 5) 8/14 (ix4 S 8ln(7max,io))~ 1/7 G. (35) 

Therefore, the discussion above reveals another important 
physical aspect of 5 m i n ; it limits the range of B- values among 
which one should choose in order to obtain the minimum 
power of the jet. If 8 — <5 m i n the dependence of Pj°t s on B 
can also be derived by combining equations (|15|l . (|24p and 

(ED 

P&(B; S min ) oc ai B- 4 + a 2 B 3 for B < B q (36) 
and 

P$ B (B;5 mia ) ochB 12 ' 7 + b 2 B- 9/7 for B > £> q , (37) 

where a\ t 2 and b\ t 2 are constants. Note that the functional 
form of -Fjet s is n °t trivial and it consists of four power-law 
segments. We will return to this point in the next section, 
where we perform the numerical approach in fitting the 7- 
ray observations of 3C 279. 



3 NUMERICAL RESULTS 

In what follows we will present a method of fitting the Febru- 
ary 26, 2006 high energy observations of quasar 3C 279 us- 
ing an one-zone hadronic model. Our working framework is 
analogous to that adopted in the previous section but with 
two main differences: 

(i) no assumptions about the relative importance of syn- 
chrotron cooling with respect to that due to photohadronic 
processes for the proton distribution are made; particle 
and photon distributions are self-consistently obtained as 
the solution to a system of three coupled integrodiffer- 
ential equations. This is done with the help of the nu- 
merical code descr i bed i n iMastichiadi s fc Kirkl l|l995l ) and 
iMastichiadis et al.l (|2005h . 

(ii) x 2 statistics was used to determine which proton- 
synchrotron spectra provide 'good' fits to the VHE 7-ray 
data. 

In order to minimize the number of free parameters we keep 
fixed the following: 7 m i n = 1, s = 2 and R — 3 x 10 16 cm. We 
note that in the analytical formulation of section 2, we have 
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used explicitly the parameter R, although it combines an ob- 
servable quantity, i.e. the variability timescale t V ar with the 
Doppler factor 8. For this, the derived values of 8 through- 
out the present work should be checked a posteriori against 
the relation imposed by the variability of the source 

s> 4R»(l + z) t (3g) 

^days 

where idays is the observed variability timescale normalized 
to 1 day. One could instead work with the observable t vaT 
and incorporate this extra constraint in his/her analytical 
treatment. Possible effects of different adopted values for R 
and s on our results will be discussed in section 4. 

3.1 The method 

Our aim is to produce a parameter space which gives 'good' 
fits to February 2006 observations of 3C 279. Application 
of existing theoretical models to AGN observations leads to 
one set of parameters that minimizes the x 2 -value. However, 
the goodness-of-fit may pose some interesting questions, es- 
pecially when two very different sets of parameters may have 
very similar values of \ 2 - For this reason, in what follows, 
we do not restrict ourselves to the best fit with the mini- 
mum x 2 -value, but we rather relax the definition of a 'good' 
fit. Thus, fits to the TeV data that have x?ed < 1-5 are 
characterized as 'good' and they are obtained for different 
combinations of the parameters. This, as we will show in the 
next section, results in the formation of a parameter space 
instead of a single set of accepted parameter values. 

Since we keep 7 m i n , s and R fixed the number of free 
parameters in the context of a pure hadronic model reduces 
to four: B, 7 max , £ P and S. Thus, we search for combinations 
of the aforementioned parameters that provide good fits to 
the TeV data. At the same time we have treated the X-ray 
observations as upper limits: As long as the proton induced 
emission is below the X-ray data, we assume that one can 
always find a fit to them by using a suitably parametrized 
leptonic component. On the other hand, if the emission due 
to photon quenching is above the X-rays, then we discard 
the fit. The steps of the algorithm followed are 

(i) Choose a value for the magnetic field strength B. 

(ii) Choose a value for the maximum proton Lorentz fac- 
tor 7max starting from the highest possible value, which is 
imposed by the Hillas criterion (see equation (|14|) '). 

(iii) Choose an injection compactness for the proton dis- 
tibution £ p . For high enough values, automatic quenching 
sets in and creates a soft-photon component that exceeds 
the X-ray observations. 

(iv) Choose a value for the Doppler factor 8 after taking 
into account that the variability timescale of 3C 279 does 
not exceed that of one day. 

For each triad (B, 7 max , l P ) two values of 8 can be found 
that correspond to fits with Xied ^ 1-5. This can be easily 
explained by the parabolic shape of the Xrcd curves. Fig- 
ure [5] shows Xrcd as a function of 8 for B = 40 G and 
7max = 4 x 10 9 . Different curves correspond to different val- 
ues of £ p . The part of the curve that lies below the horizontal 
line with Xrcd = 1-5 provides good fits. Its projection on the 
horizontal axis defines an interval of <5-values. For simplify- 
ing reasons, we consider that each monotonic branch of the 



3.5 



3 - 




1.2 1.4 1.6 1.8 2 2.2 

log 6 

Figure 2. Reduced \ 2 as a function of the Doppler factor S for 
B = 40 G, 7 ma x = 4 X 10 9 and £ p = 10" 8 (solid line), £ p = 1CT 7 
(dashed line), £ p = 10 -6 (dotted line), £ p = 6.3 X 10 -6 (dashed- 
dottcd line) and £ p = 2.5 X 10~ 5 (dashed-double dotted line). The 
horizontal solid line with Xred shows the upper limit below which 
a fit is characterized as 'good'. The line X 2 c( j = 1 is also shown. 



X 2 ed- cu rve that lies below the horizontal line is represented 
only by one 5-value, i.e. only by one point on the horizontal 
axis. Thus, for each £ p we find two representative values of 
the Doppler factor. The error we make in this case is not 
large since the curves are very steep. Note also that Doppler 
factors lying between the two representative 5-values also 
provide good fits. If we repeat the above proceedure for var- 
ious 7max we can then create a parameter space for £ p and 
8 for a specific value of the B-field. An example is shown in 
figure[3]for B = 40 G and different 7 m ax- For the case consid- 
ered here, 7 max ranges from 10 9 ' 6 (horizontally dashed area) 
to 10 10 (vertically dashed area) with a step of 0.2 in loga- 
rithmic units. The envelope of each striped area is the result 
of the two representative values of 8 defined for each £ p , as 
previously described. Note that in a log-log plot the lines 
of the envelope are power laws with an exponent —0.22, i.e. 
8 oc £p °' 22 . This almost coincides with the value —2/7 found 
in our analytical approach (see equation lf 3[) . Values of 7 max 
outside this range do not provide good fits and are therefore 
discarded. Thus, the shaded areas in Figure [3] depict all the 
possible combinations of £ p , 7 max and 8 that provide good 
fits to the TeV observations. We note that one can use low 
values of the proton injection compactness £ p and get ac- 
ceptable TeV fits but this can be done only if he/she allows 
the Doppler factor to take high values - this is due to the 
fact that the observed and blob frame luminosities are re- 
lated by Lobs oc <5 4 Z/i nt , while L- m % oc L p oc £ p . On the other 
hand, higher values of £ p naturally result in lower 8 values. 
However, as stated earlier, we find that 7— ray quenching 
plays an important role when £ p takes higher values, fn this 
case the spontaneously produced soft photons increase and 
eventually will start violating the X-ray observations, de- 
stroying the goodness of the fit. Therefore, an upper limit is 
imposed on the allowed values of £ p which implies in turn, 
the existence of a minimum value of the Doppler factor (see 
equations ((23)) and pi))). 
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Figure 3. Parameter space of pairs (£ P ,S) for B = 40 G and 
different 7 max that correspond to TeV fits with reduced \ 2 < 1-5. 
Each shaded area corresponds to a different value of 7 max that 
ranges from 10 9 ' 6 to 10 10 with a step of 0.2 in logarithmic units. 
For each value of l p two values of <5 can be determined, which form 
the envelope of each striped area. The shaded areas enclosed by 
the solid lines correspond also to values of (£ p , 8) for which food 
fits to the TeV data (Xred < 1-5) can De obtained. 



3.2 Minimum Doppler factor 

Figure [4] exemplifies the effect of 7-ray absorption on the 
adopted values of £ p and 5 for fitting the TeV data. The 
MAGIC TeV data points show the measured flux cor- 
rected for intergalactic 77 absorption. The lowest possible 
level of extragalactic background light (EBL) according to 
iPrimack etafl |2005l ) has been used for the de-absorption of 
the VHE 7-rays. The full line curve depicts the last good- 
fit spectrum we obtain for a high t p . Here t p = 4 x 10 , 
7max = 10 9 ' 8 and 5 m i n = 18.6. Other parameters used for 
this plot are: 7 min = 1, R = 3 X 10 16 cm, B = 40 G and 
s = 2. Note first that the minimum Doppler factor derived 
here does not violate the inequality (|38[) for tdays = 1, and 
second that B 2 7 max satisfies the constraint given by relation 
(|6|. The non-linear cascade due to quenching has produced 
a soft component that is almost at the observed X-ray flux. 
Any attempt to increase £ p further will produce a higher flux 
of soft photons which would violate the X-rays. To demon- 
strate this point further, we have ran the code for a high 
value of £ p = 10 -3 . This value is well inside the quenching 
regime and thus should be excluded due to strong X-ray pro- 
duction. However, by turning artificially 77 absorption off, 
we inhibit the growth of the non-linear cascade and there- 
fore a good fit can be obtained with a 5 as low as 10 (dashed 
line). Moreover, if 77 absorption is treated in an approxi- 
mate semi-analytic manner as a linear absorption process, 
automatic quenching of 7-rays will not even occur , since 
it is a purely non-linear absorption process. Therefore, no 
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Figure 4. Multiwavelength spectra of 3C 279 in the context of a 
pure hadronic model for R = 3 X 10 16 cm, B = 40 G, 7 m i n = 1, 
7max = 6.3 X 10 9 and s = 2.0. The solid curve is obtained with 
all processes taken into account for £p — 10 ^' and (5 m i n — 18.6, 
whereas the dotted one is obtained for £ p = 10~ 3 and S = 10 
after artificially ignoring 77 absorption. Squa res represent the 
VHE-detection by MAGIC llAlbert et al.ll2008T) after correcting 
for intergalactic 77 absoprtion, the bowtie represents the Swift 
data and filled circles are IR d ata during the flare of June 1991 
taken by llHartman et al,|ll996T) . 
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Figure 5. Minimum value of the Doppler factor <5 m i n (points) as 
a function of the magnetic field strength B. The dashed line is 
the result of spline interpolation between the points. 



limiting value of the proton and 7-ray compactness would 
be found and the transition of the hadronic system to su- 
percriticality would not be seen. As a result one would find 
erroneously a more extended parameter space that fits the 
data. In such case, the Doppler factor would not be lim- 
ited to a lower value apart from the one found by the usual 
gyroradius arguments made in section 2. 

As a next step we have obtained <5 m i n for various val- 
ues of the magnetic field B. This is depicted in figure [5j 
where the two branches derived analytically in the previous 
section (see figure [TJ are clearly seen. The power law depen- 
dence of (5 m i n on B can be modelled as of <5 m i n oc B~ a with 
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a ~ 2.85 for the low-B branch and a ~ 0.05 for the high-B 
branch. We note that the numerically derived power-law ex- 
ponents of the two branches are very close to those given by 
equations (|15[) and (|24[) respectively. A new feature of figure 
[5] is that for very strong magnetic fields, where synchrotron 
cooling affects a significant part of the proton power-law dis- 
tribution, the minimum Doppler factor required for a good 
fit slightly increases with increasing B. 



3.3 Energetics 

The previous analysis makes evident that acceptable fits to 
the TeV observations of 3C 279, within the context of the 
hadronic model, can be found for a large choice of mag- 
netic fields, maximum proton energies and injection lumi- 
nosities. We note also that these values are similar to the 
ones found in standa rd hadronic modelling of TeV sources 
ijBottcher et al.ll2009l l. 

Having already derived <5 m i n for each magnetic field 
strength, we numerically then derive the steady state proton 
energy density m* = M p (c5mm(-B)) and plot this quantity as 
a function of B - see figure [6] The system becomes magnet- 
ically dominated either for high enough (B > 10 G) or low 
enough (B < 0.7 G) magnetic fields. These values should be 
compared to Bh and B\ given by equations (|34|l and (|35[) 
respectively. The two curves intersect at two points, exactly 
as was analytically derived (see section 2.2). At these points 
the two energy densities are equal and this corresponds to a 
minimum of the total energy density, i.e. to an equipartition 
between the particles and the magnetic field - we empha- 
size that the quantity Up is not the minimum particle energy 
density. The fact that <5 m i n has a strong dependence on B 
at low values is also reflected to lit. Our analytical result 
given by equation ()28|) should be compared to the numeri- 
cal one, i.e. up~ tx B 8 - 7 . For B > 2 G (c.f. B q of equation |27|) 
the particle energy density decreases with increasing mag- 
netic field as B~ 14 . This trend was also found analytically. 
The power-law exponent however is slightly different - see 
equation l(2"9")) . 

From equation ((30)) and for 8 = 8 m i n (B) and £ p = 
^p(<5 m in) we have calculated numerically the observed jet 
power, which becomes a function only of the magnetic field, 
i.e. P°^ s ((5 m i I1 (B)). This is shown in figure [3 The function 
Pj°t s (<5 m in(B)) shows two local minima for two values of the 
magnetic field that differ more than one order of magni- 
tude. For a wide range of B values the calculated jet power 
is rather high ~ 10 48 — 10 49 erg/s in comparison to lep- 
tonic models (|Celotti fc Ghisellini 2008). In a log-log plot 
the curve consists of four distinct power-law segments la- 
beled with the numbers 1 to 4. In Table 1 the power-law 
exponents of the four segments and the corresponding ana- 
lytically derived values given by equations ([36j) and (|37[) are 
listed. Apart from the first segment both results are in good 
agreement. 

We have already shown that the jet power minimizes 
whenever 8 and B satisfy equation ()32[) . In this case the 
system is close to a state of equipartition. We proceed to 
investigate whether these results are verified numerically. 
For each value of B, we therefore search among the triads 
(7max, 8) found during the fitting process for that specific 
set of values that minimizes the jet power. Figures [8] and [9] 
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Figure 6. Logarithmic plot of the proton energy density in the 
case of the mimimum Doppler factor u* (points) and of the mag- 
netic energy density (solid line) as a function of the magnetic field. 
Spline interpolation between the points results in the dashed line. 




Figure 7. Jet power as a function of the magnetic field for <5 = 
5 m i n (points). The dashed line is the result of spline interpolation. 
Numbers label the different power-law segments of the curve. 



exponent 


Number of segment 


numerical 


analytical 


(1) 


-2.50 


-4.00 


(2) 


+2.30 


+3.00 


(3) 


-1.30 


-1.28 


(4) 


+1.70 


+1.71 



Table 1. Power-law exponent of P°^ s (5 m m(B)) as derived from 
the numerical and analytical treatment. 
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Figure 8. Minimum jet power value (open circles) as a function 
of the magnetic field strength. The horizontal line denotes the 
analytically derived minimum value of the jet power (equation 
I33P and the dashed line, which is same as in figure [jj is shown 
for comparison reasons. Open circles that lie on the dashed line 
correspond to the jet power calculated for <5 m ; n , since <5 < 5 m i n 
would be required to minimize the jet power. 



Figure 9. Proton energy density (open circles) as a function 
of the magnetic field strength, for those parameter sets that cor- 
respond to the minimum jet power. The magnetic field energy 
density is shown with the solid line. The remaining symbols are 
the same as in figure [8] 



and the corresponding proton energy density as a function 
of B. Dashed lines are same as in figures [6] and [7] and are 
shown for comparison reasons. 

Figure [8] shows that for 0.7 < B sC 10 G all points 
lie close to the minimum value -Pj°tf m i n — 6.5 x 10 47 erg/s, 
whereas for B > 10 G and B < 0.7 G the minimum jet power 
found from our numerical data sets coincides with that cal- 
culated for 8 — <5 m i n (open circles lying on the dashed line). 
As already discussed in section 2.3 the jet power calculated 
in this case is not the lowest possible minimum value, since 
this would be obtained for Doppler factors less than 5 m in- 
Note that the transitions occur at B « B\ w 0.7 G and 
B « B h w 10 G. 

In figure [9] one can see that u p f%s ub for Bi ^ B ^ Bh, 
which is in complete agreement with the analytical estimates 
of section 2.3 . The fact that the calculated jet power is not 
the minimum possible one, can also be verified by the fact 
that the system for B > 10 G and B < 0.7 G is far from an 
equipartition state. 

Each of the points shown as open circles in figures [8] and 
[9] corresponds to a set of parameters (B, 7 max , i P , 8) with 
8 ^ 5 m in, for which the jet power is actually the minimum 
possible in the context of a hadronic model. In figure \W\ 
we plot the ratio <5/<5 m i n against the magnetic field strength. 
For the range 0.7 < B ^ 10 G where the system is close to 
equipartition and the power of the jet is minimum, the ratio 
becomes larger than unity. 

The conclusions drawn from the figures of the present 
section can be summarized as follows: if one chooses to 
model the VHE 7-ray spectra using values of (B, 7 max , £ P , 8) 
that minimize both the required jet power and the total 
energy density of the system, i.e. model the system using 
optimum energetic conditions, then one will be confronted 
with the requirement of a high Doppler factor. On the other 
hand, acceptable fits of the TeV data using the minimum 
possible Doppler factor lead to a particle dominated system 
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Figure 10. Ratio 5/i5 m i n as a function of the magnetic field. Here 
8 is the Doppler factor required for minimization of the jet power. 
Interpolation between the points results in the dashed line. 



and to high values of P£, t B for a wide range of B values, i.e. 
the energetic requirements in this case are higher. 



4 EFFECTS OF OTHER PARAMETERS 

In our numerical treatment we have kept fixed to specific 
values the power-law index of the injected proton distribu- 
tion, the radius of the emitting blob as well as the energy 
of the 7-ray photons e° bs . Here we discuss possible effects 
of different adopted values on our results, since our initial 
assumption could be critical for exploring the whole physi- 
cally allowed parameter space for the system of 3C 279 at 
the observed state. 
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4.1 Power-law index s 

Throughout the present work we have assumed that the pro- 
ton distribution is a power law of exponent s = 2. In this 
case the energy per logarithmic interval is the same. A dif- 
ferent adopted value would affect the proton energy density 
quantitavely and therefore our energetics estimates. Let us 
for example consider a steep power law distribution with 
s > 2. In this case more energy is injected to the protons 
with the minimum Lorentz factor, whereas protons at the 
upper cutoff of the distribution, which are responsible for the 
emitted radiation, carry only a small fraction of the energy. 
Thus, for the same proton injection compactness the total 
number N p and consequently the energy density of protons 
Up increases with increasing s. In addition, the calculated 
7-ray compactness for the same £ p is lower. Thus, a larger 
proton compactness i v and therefore energy density would 
be required in order to fit the TeV observations with the 
same Doppler factor. As an extreme example, we considered 
a very steep proton distribution with the highest value for 
the power law exponent predicted by acceleration theory, i.e. 
s — 2.5. The fit to the TeV observations was obtained using: 
B = 40 G, 8 w 25 and l v = 10~ 13 . The proton energy den- 
sity in this case is u p = 10 4 erg/cm -3 , which is much higher 
than the values presented in section 3.3. Therefore, the case 
of a fiat proton distribution is rather conservative as far as 
the energetics is concerned. 

4.2 Radius of the emitting region 

The fact that there is a very good agreement between the 
numerical results and our analytical expressions, where the 
dependence on R is explicit, makes possible the prediction 
of the effects of a different adopted value on our results. Let 
us assume a more compact source with R = 3 x 10 15 cm. 
The quantities that are directly affected by a change in the 
source size are listed below in descending order in terms of 
their dependence on R: 

(i) The minimum Doppler factor for B < B q (see 
equation I15[) . which would be increased by two orders of 
magnitude. Therefore, the 'steep' branch of the plot in 
figure [5] would be shifted upwards by a factor of two in log- 
arithmic scale. Note, that the constraint set by variability 
arguments (relation I38p will not be violated since it implies 
an even lower limit to the Doppler factor than previously. 

(ii) The equipartition magnetic field B cq given by equa- 
tion Q27p . which for a certain 8 will be increased by a factor 
of 7. 



(see equation I25[) ; the dependence on R comes through the 
constant A. Taking into account points (i) and (iv) above, 
we find that u p (<5 m i n ) oc R 4 for B < B q and u p (5 m i n ) tx R~ 4 
for B > _B q . Since B q also increases by a factor of 4, one 
expects to find significantly smaller values for M p (<5 m in) than 
those shown in figure O for a wider range of B-values (up 
to ~ 14 G). The jet power obtained in this case has also 
a strong dependence on R, which is embedded in the con- 
stants ai,2 and bi^ of equations (|36|) and (|37[) respectively. 
More specifically one finds that ct\ oc R~ 2 and «2 oc R 2 
whereas b\ oc R 10 ^ 1 and hi oc i? _4//r . The shape of the curve 
shown in figure [7] would be transformed by shifting the dif- 
ferent power-law segments vertically and horizontally, since 
the values of the magnetic field where the local minima oc- 
cur would also be affected (see equations [M] and I35p . It is 
important, however, to note that the minimum Doppler fac- 
tor in this case would be extremely high (see also point (i) 
above). This makes the scenario of a more compact 7-ray 
emitting region for 3C 279 less plausible. 



4.3 Energy of 7-ray photons 

Another question that naturally arises is how our results 
would change if the fitting method was used for another set 
of 7- ray observations, e.g. in the GeV regime. Note that 
some rec ent contemporary 7- and X-ray observat ions of 3C 
279 (e.g. lAbdo et al.ll2010l ; iHavashida et alj|2012h could be 
an interesting case. In principle, the method outlined in sec- 
tions 2 and 3 can be applied with the only difference that the 
effects of quenching will not be seen, as automatic quench- 
ing cannot set in for lower 7-ray energies at least for typical 
magnetic field strengths. This is examplified in figure 111! 
where characteristic values of the magnetic field are plot- 
ted against 7 max in the extreme case of 8 = 1. Specifically, 
the thick solid line corresponds to B q (see equation (|19[) ) 
and divides the parameter space into two regions. For B 
values that lie above this line, the feedback criterion is sat- 
isfied. The thin solid line corresponds to the Hillas criterion, 
i.e. B ^ (m p c 2 /e_R)7 max . Therefore, the parameter space 
below this line is not allowed. Finally, the thin and thick 
dotted lines show the locus of B and 7 max values, that cor- 
respond to e° bs equal to 0.1 TeV and 1 GeV respectively, 
i.e. B — e° bs (m p c/fie)7max<5 -1 for 5 = 1. Therefore, in the 
case where GeV observations were used, one would need ex- 
tremely high values of the magnetic field, in order to see the 
effects of automatic quenching. These would be even higher 
if one would take into account the exact value of the Doppler 
factor. 



(iii) The magnetic field above which the feedback crite- 
rion is satisfied. We remind that the scaling of B q oc R~ 3 ^ 5 
has been derived after taking into account the constraint of 
the Hillas criterion on y max . 

(iv) The minimum Doppler factor set by the automatic 
quenching is rather insensitive to changes of R and it would 
be increased just by a factor of 2 (see equation Ij24| l), 

The proton energy density as well as the jet power calcu- 
lated in the case of a blob moving with S — Smin will be also 
affected by a change in R. We remind that tt p oc R~ 3 8~ 7 ^ 2 



5 DISCUSSION 

Hadronic models have been used extensively for fitting the 
7-ray emission from Active Galactic Nuclei. Usually, detailed 
fitting to the multiwavelength spectrum of these objects re- 
quires, in addition to a population of relativistic protons, 
the presence of an extra leptonic component which is re- 
sponsible for the emission at lower energies (radio to UV or 
X-rays). 

As it was shown recently ( Stawarz fc Kirkl 120071 ; 
iPetropoulou fc Mastichiadis! l201ll . I2OI2I ) compact 7-ray 
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Figure 11. Log- log plot of characteristic values of the magnetic 
field against the maximum proton Lorentz factor. The space be- 
low the thin solid line is not allowed, since it violates the Hillas 
criterion. The feedback criterions is satisfied for values of B and 
7max that lie above the thick solid line, which represents B q . Fi- 
nally, the thin and thick dotted lines B and 7 max values, that 
correspond to e° bs equal to 0.1 TeV and 1 GeV respectively. 



sources can be subject to photon quenching which can re- 
sult, if certain conditions are met, to automatic 7-ray ab- 
sorption and redistribution of the absorbed 7-ray luminosity 
to electron-positron pairs and lower energy radiation. This 
is a non-linear loop which can operate even in the absence 
of an initial soft photon population and is expected to have 
direct consequences on the aforementioned models. For the 
present application other feed back loops, such as the Pair - 
Production-Synchrotron one l|Kirk fc Mastichiadid 1 19921) , 
are less relevant as they usually operat e at lower proton en- 
ergie s and at higher energy densities (|Dimitrakoudis et al.l 
l2012h . 

Aim of the present paper is to set a general framework 
for investigating the effects of photon quenching on the pa- 
rameter space available for modelling 7-rays in the context 
of a hadronic model. As an example we concentrated on the 
February 2006 observations of the blazar 3C 279 which were 
performed simultaneously in the TeV and X-ray regimes. 
To this end we have focused only on fitting the high energy 
spectrum using a synchrotron proton emission and treating 
conservatively the X-ray observations as an upper limit. We 
found that for a wide range of parameters, automatic pho- 
ton quenching plays a crucial role as its onset produces X- 
rays which violate the observations. To be able to assess its 
impact on the parameter space, we have relaxed the usual 
goodness-of-fit method by accepting fits with Xicd < 1-5. 
Our results indicate, in agreement with other researchers in 
the field, that the hadronic model requires in general high 
magnetic fields. Additionally we find that the presence of au- 
tomatic quenching limits the proton luminosity (or, equiva- 
lently the proton energy density in the case of non substan- 
tial proton cooling) which, in turn, results in a minimum 
value of the Doppler factor <5 m in- Interestingly enough, these 
ideas do not apply to leptonic models because they favour 



much lower values for the magnetic field. For these values 
photon quenching does not operate, since its feedback cri- 
terion is not satisfied. The latter is a necessary condition 
for the emergence of automatic quenching and requires for a 
given 7-ray energy a certain value of the magnetic field B q , 
for the absoprtion loop to operate. We have shown both ana- 
lytically and numerically, (see figures [1] and [5] respectively), 
that the minimum Doppler depends on the magnetic field 
strength in a different way, depending on the relative rela- 
tion of B and B q . Specifically, if the magnetic field strength 
is above S q , then <5 m i n oc _B _1//r . For values of the mag- 
netic field with B < B q quenching is not relevant. However, 
we have shown using arguments based on the particles gy- 
roradii, that also in this case a lower limit to the Doppler 
factor exists, which depends strongly on B, i.e. 8 m in oc B~ 3 . 
Therefore if one wants to adopt a low magnetic field for the 
radiating region, fits to the TeV 7-rays are still possible but 
at the expense of a very high value of the Doppler factor. 

The fact that quenching does not allow the Doppler fac- 
tor to become smaller than some value is intriguing and leads 
naturally to the investigation of the proton energy density 
inside a blob which moves with this characteristic value. 
In this case, we showed that there are two values of the 
magnetic field that minimize the energy content, one corre- 
sponding to the S branch with B < B q and the other to the 
one with B > _B q . For values of the magnetic field between 
the two equipartition magnetic fields, the emission region is 
particle dominated and the ratio u p /ub can be as large as 
10 3 (see figure Note however, that if a Doppler factor 
twice the minimum one is adopted, then the calculated pro- 
ton energy density for the same magnetic field will be lower 
by almost an order of magnitude - see equation (I12|l . Fur- 
thermore, we have calculated the jet power in the case where 
S — <5min and shown both analytically and numerically that 
it is a function of B. The jet power is rather high (10 47 -10 49 
erg/sec) for the whole range of B values, as expected in the 
context of hadronic modelling. 

We have repeated the above calculations by relaxing 
the condition 8 = S m i n while requiring the parameters to 
be such as to minimize the power of the jet. In this case 
we have shown that the energy content of the blob is also 
minimized (see figures [8] and [9} . For adopted B values lying 
between the two equipartition values, the jet power can be 
minimized, albeit at the cost of a high 8 value. However, 
we have shown that for magnetic field values outside of this 
range, jet power minimization is not possible as this can be 
achieved only if 8 < 5 m i n - Thus, the existence of a minimum 
value for 8 has indirect implications on the energetics of the 
system. 

An interesting question is whether a detailed fitting to 
the X-ray observations of 3C 279 with the addition of an 
extra leptonic component would bring any change to the 
basic ideas presented here. In this case apart from auto- 
matic quenching, the linear absorption of 7-rays on the X- 
ray photons emitted by the leptonic component, is also at 
work. However, by including this component and repeating 
our numerical calculations of §3, we found that our results do 
not change. This is due to the fact that the X-ray luminosity 
of 3C 279 is rather low and therefore the effects of linear 77 
absorption are minimal, at least up to the compactnesses 
above which the automatic photon quenching sets in. Apart 
from the absorption of 7-rays on the synchrotron photons 
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emitted by the 'extra' leptonic component discussed above, 
inverse compton scattering of these photons to higher ener- 
gies by the same leptonic component would be an additional 
mechanism at work. The upscattered photons would lie in 
the hard X-ray and 7-ray energy range and would affect 
our calculations only if their luminosity L ssc would be com- 
parable to that carried by the synchrotron component L syYL . 
While estimating the ratio u syn /uB for a wide range of B and 
<5mm values used in the present work (see figure [SJl we find 
that this is always much smaller than unity, and therefore 
the inverse Compton scattering of an extra leptonic compo- 
nent which would fit the radio to X-ray observations does 
not interfere to the quenching mechanism. 

The estimation method proposed in the present work 
can be regarded as an extension of the widely used method 
for estimating the equipartition magnetic field using radio 
observations. In our case, the leptonic synchrotron compo- 
nent is replaced by the proton synchrotron emission and the 
radio by the VHE 7-ray observations. The innovative feature 
of our method is the estimation of a minimum Doppler fac- 
tor, which is the result of automatic photon quenching. This 
can be more of relevance to TeV observations rather than the 
Fermi ones as quenching at GeV energies requires very high 
values of the magnetic field - see figure [TT] The fact that 
our numerical results are in very good agreement with the 
analytical calculations offers a fast yet robust way for esti- 
mating many physical quantities of the 7-ray emitting region 
of TeV blazar jets. We caution however the reader, that the 
effects of automatic quenching, i.e. when no soft photons are 
initially present in the emitting region, can be seen only in a 
self-consistent treatment of the radiative transfer problem. 



Hillas A. M., 1984, ARAA, 22, 425 
Kirk J. G., Mastichiadis A., 1992, Nature, 360, 135 
Mannheim K., Biermann P. L., 1992, A&A, 253, L21 
Maraschi L., Ghisellini G., Celotti A., 1992, ApJ, 397, L5 
Mastichiadis A., Kirk J. G., 1995, A&A, 295, 613 
Mastichiadis A., Protheroe R. J., Kirk J. G., 2005, A&A, 
433, 765 

Miicke A., Protheroe R. J., 2001, Astroparticle Physics, 15, 
121 

Miicke A., Protheroe R. J., Engel R., Rachen J. P., Stanev 

T., 2003, Astroparticle Physics, 18, 593 
Murase K., Dermer C. D., Takami H., Migliori G., 2012, 

ApJ, 749, 63 
Pacholczyk A. G., 1970, Radio Astrophysics 
Petropoulou M., Mastichiadis A., 2011, A&A, 532, A11+ 
Petropoulou M., Mastichiadis A., 2012, MNRAS, 421, 2325 
Primack J. R., Bullock J. S., Somerville R. S., 2005, in 

American Institute of Physics Conference Series, Vol. 745, 

High Energy Gamma-Ray Astronomy, Aharonian F. A., 

Volk H. J., Horns D., eds., pp. 23-33 
Protheroe R. J., Miicke A., 2001, in Astronomical Society 

of the Pacific Conference Series, Vol. 250, Particles and 

Fields in Radio Galaxies Conference, Laing R. A., Blun- 

dell K. M., eds., p. 113 
Stawarz L., Kirk J. G., 2007, ApJ, 661, L17 
Tavecchio F., Ghisellini C, Bonnoli G., Foschini L., 2011, 

MNRAS, 414, 3566 



6 AKNOWLEDGEMENTS 

This research has been co-financed by the European Union 
(European Social Fund - ESF) and Greek national funds 
through the Operational Program "Education and Lifelong 
Learning" of the National Strategic Reference Framework 
(NSRF) - Research Funding Program: Heracleitus II. In- 
vesting in knowledge society through the European Social 
Fund. We would like to thank Dr. Anita Reimer and the 
anonymous referee for usefull comments on the manuscript. 



REFERENCES 

Abdo A. A. et al., 2010, Nature, 463, 919 
Albert J. et al., 2008, Science, 320, 1752 
Bottcher M., Reimer A., Marscher A. P., 2009, ApJ, 703, 
1168 

Celotti A., Ghisellini G., 2008, MNRAS, 385, 283 
Chatterjee R. et al., 2008, ApJ, 689, 79 
Dermer C. D., Schlickeiser R., 1993, ApJ, 416, 458 
Dermer CD., Schlickeiser R., Mastichiadis A., 1992, A&A, 
256, L27 

Dimitrakoudis S., Petropoulou M., Mastichiadis A., 2012, 
International Journal of Modern Physics Conference Se- 
ries, 8, 19 

Harris D. E., Krawczynski H., 2002, ApJ, 565, 244 
Hartman R. C. et al., 1996, ApJ, 461, 698 
Hayashida M. et al., 2012, ArXiv e-prints 



© 2012 RAS, MNRAS 000.1111121 



